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We discuss the amplitudes describing A^-gluon scattering in type I superstring theory, on a disk 
world-sheet. After reviewing the general structure of amplitudes and the complications created by 



. the presence of a large number of vertices at the boundary, we focus on the most promising case of 

o , 

\^ , maximally helicity violating (MHV) configurations because in this case, the zero Regge slope limit 
O ' 

(a' — > 0) is particularly simple. We obtain the full-fledged MHV disk amplitudes for = 4, 5 and 

^ • N = Q gluons, expressed in terms of one, two and six functions of kinematic invariants, respectively. 

r~| , These functions represent certain boundary integrals - generalized Euler integrals - which for A^ > 6 

^ correspond to multiple hypergeometric series (generalized Kampe de Feriet functions). Their a'- 

^ • expansions lead to Euler-Zagier sums. For arbitrary A^, we show that the leading string corrections 



to the Yang-Mills amplitude, of order ©(a'^), originate from the well-known a'^Tri^^ effective 
interactions of four gauge field strength tensors. By using iteration based on the soft gluon limit, 
we derive a simple formula valid to that order for arbitrary A^. We argue that such a procedure can 
be extended to all orders in a' . If nature gracefully picked a sufficiently low string mass scale, our 
results would be important for studying string effects in multi-jet production at the Large Hadron 
Collider (LHC). 

PACS numbers: 11.25.Wx, 11.30.Pb, 12.38.Bx 
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I. INTRODUCTION AND REVIEW 



Since Rutherford's times, elementary particle physics relies on scattering experiments. 
The physical cross sections, determined by the scattering amplitudes, reflect the properties 
of underlying interactions. In the framework of the standard model, high energy scatter- 
ing experiments allow probing inside hadrons, into the gauge interactions of quarks and 
gluons. Due to the asymptotic freedom of Quantum Chromodynamics (QCD), the corre- 
sponding amplitudes can be computed perturbatively, order by order in the QCD coupling 
constant. Already at the tree level, such computations can be quite complicated, especially 
when a large number of external particles is involved, like in the scattering processes de- 
scribing multi-jet production at hadron colliders. After more than thirty years of steady 
progress in perturbative QCD, we have a very good understanding of the tree- level scatter- 
ing, completely sufficient for studying jet physics in the upcoming Large Hadron Collider 
(LHC) experiments at CERN. Hopefully, LHC will reach beyond the standard model, and 
the signals of new physics will rise above the QCD background. 

Although many scenarios have been proposed beyond the standard model, there is no 
clear prediction for the energy scale of new physics. Even if no spectacular effect like, say, a 
direct production of Kaluza-Klein particles, is discovered at LHC, some sub-threshold effects 
could be observed, due to the presence of contact interactions induced by virtual particles 
too heavy to be produced on-shell. In this paper, wc investigate such effects in superstring- 
based scenarios, where the scale of new physics is determined by the Regge slope a' of 
mass dimension —2. Massless gauge bosons are separated by a mass gap of l/yfa! from the 
massive string modes. Traditionally, the Regge slope and the respective string mass scale 
had been tied to the Planck mass, however more recently, some serious consideration was 
given to D-brane models with much lower string mass scale, possibly even within the reach 
of LHC [1, 2]. The full-fledged string amplitudes depend on resulting in large corrections 
to Yang-Mills (YM) amplitudes once some some kinematic invariants characterizing energy 
scales involved in the scattering process become comparable to 1/ -\/c7. 

We work in the framework of open type I superstring theory compactifled to four dimen- 
sions, with gluons being open string excitations. In the tree approximation, the multi-gluon 
amplitudes are computed on a disk worldsheet, with the vertices inserted at the boundary. 
They do not depend on the compactiflcation manifold because they are completely de- 
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termined by two-dimensional superconformal field theory describing four space-time string 
coordinates of the worldsheet. In particular, it does not matter whether supersymmetry is 
broken or not by compactification. Note that the a' — > limit is completely determined by 
pure Yang-Mills theory. 

An important feature of open string (disk) computations is that they yield gluon ampli- 
tudes in a very particular, color-decomposed form: 

A'''\{ki,X^,a^})=g''-^ J2 Tr(r»<'W---r«-w) ^(a(l^^),...,a(7V^^)) , (1) 

2 

where g is the gauge coupling (|^ = ct^), ki,Xi are the gluon momenta and helicities, 
and T"' are matrices in the fundamental representation of the gauge group [3], describ- 
ing the color states of N gluons. We consider amplitudes with all momenta directed in- 
ward. Each color trace factor is associated by Chan-Paton rules to one partial amplitude 
A{a{l^^), . . . , (t{N^'^)) containing all the kinematical information. Sn is the set of all per- 
mutations of objects, while Zjv is the subset of cyclic permutations, which preserves the 
trace; one sums over the set Sf^/'L^ in order to include all orderings of gluon vertices, sweep- 
ing out all distinct cyclic orderings in the trace. A similar color decomposition is routinely 
used for multi-gluon amplitudes in QCD [4, 5]. 

In QCD, there exists a subclass of amplitudes that arc described, at the trcc-lcvcl, by a 
simple analytic formula valid for arbitrary number A^ of gluons. Assume that two gluons, 
with the momenta ki and A;2, in the color states described by the matrices and T"^, 
respectively, carry negative helicities while the rest of gluons, with the momenta and color 
charges (/cs, T"-^), . . . , [kN, T"^), respectively, carry positive helicities. Then the partial am- 
phtude for such a "maximally hehcity violating" (MHV) configuration, associated to the 
Tr ( • • • T"^) Chan-Paton factor, is given by [6, 7] 

where {ij) are the standard spinor products associated to the momenta ki, kj, in the notation 
of [4, 5]. Other partial amplitudes can be obtained from Eq.(2) by applying appropriate 
permutations to the cyclic denominator (12)(23)(34) ■ ■ ■ (A^l). For example. 
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Thus the full MHV amplitude is 

2-. 3^ . . . . N^) ^ .."-(12)-^^E^ ^^^,^^7,^'3^;;;,J°:>^^^^^ . (4) 

where v = cr{i). The origin of the striking simphcity of MHV amphtudes is not clear. 
Most hkely, it is related to some (partial) integrability properties of QCD. There is also an 
interesting duahty between Yang-Mills theory and twistor strings [8] that led to a new inter- 
pretation of Eq.(2) and turned out very useful for developing more efficient computational 
techniques in perturbative QCD [9, 10]. The amplitudes describing non-MHV helicity con- 
figurations are known to be more complicated. Thus when studying multi-gluon scattering 
in string theory, it is natural to use MHV configurations as a starting point, in order to 
find out if the full-fledged string amplitudes can also be described by some simple analytic 
formulas valid to all orders in a'. Our results show that this is indeed the case. 

In order to write down the amplitudes in a concise way, it is convenient to introduce the 
following notation for the kinematic invariants characterizing A/^-particle scattering: 

bin = i^i + ^i+l + ■■■ + ki+n? , (5) 

e{i,j,m,n) = a''^ eaPf,„ k°'k'^kX (6) 

where ki denotes the momentum of i-th particle, with the cyclic identification i + N = i, 
and Capij,!^ is the four-dimensional Levi-Civita symbol. The factors of a' render the above 
invariants dimensionless. Note that the momenta are always subject to the momentum 

N 

conservation constraint, ^ — 0, and all gluons are on-shell, kf — [ij^ = 0. It is also 

i=l 

convenient to introduce 

Sij = 2a' kikj . (7) 

By using momentum conservation, these scalar products can be always expressed in terms 
of N{N — 3)/2 invariants (5). This is done for A?" < 6 in Appendix A. Note however, that 
for A^ > 6, the number of independent invariants of type (5) is smaller than A^(A^ — 3)/2, as 
explained in Section V. 

The amplitude for four-gluon scattering has been known for a long time [11-13]. All 
string effects are summarized in one Beta function (Veneziano amplitude) 

V(^\k,, k,, ks, k,) = Vi^\s,, s,) = ^^i-^'l^ ^^y.''^ , (8) 

i (^i + Si -|- S2) 



where si = = Si2, S2 = |2]^ = S23, as the formfactor of Yang-Mills amplitude: 

^(l-,2-,3+,4+) = l-(^)(si,.2)9Jt?i, (9) 

with defined in Eq.(2). An obvious but very important property of the Veneziano 

formfactor V^'^^ki, k2, k^, /C4) is its invariance under the cyclic permutations of the momenta. 
All other partial amplitudes can be obtained by applying the coset permutations a to the 
right hand side of Eq.(9), c.f. Eq.(4), now including also the cyclic formfactor. Thus the full 
four-point amplitude is 

In fact, for four (and five) gluons, all non-MHV amphtudes are vanishing [6, 14], therefore 
the above expression captures the full amplitude. It can be expanded in powers of a' by 
using 

V^^\si, S2) « 1 - y s,S2 + C(3) S1S2 (si + S2) + 0{a'^) . (11) 

The leading string correction to the Yang-Mills amplitude, which originates from the second 
term in the above expansion, of order C(q;'^), has been extensively discussed in the literature 
[11-13]. It is due to the following contact interaction term of four gauge field strength tensors: 

where the color trace is taken with the tensors F^j, in the fundamental representation. This 
interaction term will play an important role in the further discussion of A^-point amplitudes. 
Formally, it can be obtained from the 0{a'^F^) term appearing in the low-energy expansion 
of the Born-Infeld Lagrangian of non-linear electrodynamics, by applying to it Tseytlin's 
"symmetrized trace" prescription [15]. Note that for Abelian gauge bosons, the pure Yang- 
Mills part of the amplitude (10) cancels after summing over all coset permutations and the 
low-energy expansion begins with the Born-Infeld contribution. 

The paper is organized as follows. In Section II, we give a brief description of the formal- 
ism used for calculating multi-gluon amplitudes on a disk world-sheet. Integrations over the 
vertex positions yield certain generalized hypergeometric functions, their number increasing 
dramatically with the number of gluons, therefore in addition to handling a cumbersome 



algebra, one has to figure out how to construct a proper basis of the boundary integrals. 
In Section III, we rewrite the five-gluon amplitude in an MHV form similar to four gluons, 
c.f. Eq.(9), in terms of two independent (hypergeometric) functions of kinematic invariants. 
One of them plays the role of the Veneziano formfactor, while the second is associated to 
the Xp4 contact term (12). In Section IV, we extract the MHV part of the six-gluon ampli- 
tude. Here, all kinematic information is contained in six "triple" generalized hypergeometric 
functions. We discuss the low-energy behavior of the amphtude and check that it satisfies 
all constraints based on permutation symmetries and soft/coUinear limits. We show that it 
is possible to reconstruct the result, obtained from tedious calculations, by imposing these 
constraints on a general ansatz. In Section V, we proceed to the general A^-gluon case. We 
show that all 0(q!'^) order string corrections to MHV Yang-Mills amplitudes originate from 
the interactions associated to the Xp4 effective action term. By iteration, we obtain a simple 
A^-gluon formula valid to that order and outline a recursive procedure that could make possi- 
ble a complete determination of all MHV amphtudes, to all orders in a'. In Conclusions, we 
discuss our results in a broader context of QCD and superstring theory. The paper contains 
three appendices. In Appendix A we summarize some aspects of the scattering kinematics 
for = 4, 5 and 6 gluons. In Appendices B and C we discuss various properties of the 
generalized (triple) hypergeometric functions describing the six-gluon amplitude. 
Some results of this work have been already reported in our Letter [16]. 

II. MULTI-GLUON SCATTERING ON THE DISK 

In this Section, we review the general structure of multi-gluon string amplitudes, focusing 
on the computational problems related to a large number of vertices at the boundary. It 
can be skipped by readers who are not interested in technical details. 

We are interested in superstring theory with gluons coming from open strings. A variety 
of four-dimensional models can be constructed, each of them described by a two-dimensional 
superconformal field theory (SCFT). At the disk (tree) level, the details of the "internal" part 
of SCFT associated to the compactification space do not affect the scattering amplitudes 
of four-dimensional gauge bosons. Furthermore, the entire disk boundary is attached to a 
single stack of D-branes. Thus without losing generality, we can consider type I theory with 
D9-branes and sixteen supercharges. Nevertheless, our discussion holds for both type I or 
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type IIA/B theories with Dp-branes and any gauge group. Spacetime supersymmetry can 
be preserved or broken by the internal space or by D-brane configurations. 

Gluons originate from the excitations of string space-time coordinates X'^ and their SOFT 
partners satisfying Neumann boundary conditions on the world-sheet. From all other 
SOFT fields, only the reparametrization ghost c and the scalar (j) bosonizing the superghost 
system will enter explicitly into our computations. In the (— l)-ghost picture, the vertex 
operator for a gluon with momentum A;, polarization ^ (or helicity A) and color state a is 
given by 

V^-^\z, {fc, a}) = T« e-'t'^'^ i^^iz) e'^"'"^'^ , (13) 

where z is the vertex position at the disk boundary. Note that the color state a is represented 
by the matrix in the fundamental representation of the gauge group. In the zero-ghost 
picture, this vertex operator is given by: 

l^(°)(^,{A;,e,a}) = [ dX^'{z)+i (k^) i^^iz) ] e''""^"^'^ . (14) 

The disk may be conformally mapped to the upper half plane Imz > with the real axis 
as its boundary. Hence all vertex positions are located on the real axis. The A'"-gluon disk 
amplitude is: 

/OO poo POO 

dz,^ / dZ2^ ... dZN^ 

X {V^-'\z,) V^-'\z2) V^'\z,) . . .V^'\zr,)) , (15) 

where the color part of the vertices has been factored out by following the Chan-Paton rule. 
In the above expression, Vckg is the volume of the conformal Killing group PSL{2, M) which 
leaves the boundary [Im(^) = 0] of the disk fixed. It will be canceled by fixing three positions 
and introducing the respective c-ghost correlator. Note that two vertices are inserted in the 
(—1) -ghost picture in order to cancel the background ghost charge. 

By comparing Eq.(15) with the color-decomposed form of see Eq.(l), we see that 

the partial amplitude A{a{l'^^), . . . , a{N^'^)) is obtained by integrating the correlator of the 
vertex operators over the region { — oo < Zi^ < < ■ ■ ■ < z^^ < oo }. In the following, 
we shall concentrate on the Chan-Paton factor Tr(T"i • • -T"^), i.e. in Eq.(15) we pick up 
the integration region TZ = {—oo < Zi < Z2 < . . . < zj^ < oo } and compute 

^(l^\...,7V^-) = l-CKG / {\[ dzr) (y^-'\z,) V^-'\z2) V^'\z,) . . .V^'\zm)) . (16) 
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Due to the PSL{2, IR) invariance on the disk, we can fix three positions of the vertex oper- 
ators. A convenient choice is 



Zi — —Zoo — —OO , Z2 — , ^3 — 1 , 



(17) 



which imphes the ghost factor {c{zi)c{z2)c{z3)) = —z"^. The remaining A^— 3 vertex positions 
Z4, . . . ,zn take arbitrary values inside the integration domain TZ. It is convenient to use the 
following parameterization: 

N-3 



2^4 = ^ , 2^5 = {X1X2) ^ , Zq^ ixiX2X3) ^ , 



Zn 



= n ' (18) 



Ar-3 



l+r-N 



with < < 1. The corresponding Jacobian is \dzi/dx.j\ = H ^ 

r=l 

The correlator of vertex operators in Eq.(15) is evaluated by performing all possible Wick 
contractions. It decomposes into products of two-point functions, introducing kinematic 
factors consisting of the scalar products of momentum and polarization vectors, of the form 
kikj, ^ikj and ^i^j. Schematically, one obtains 

,7 



n„ 



n 



ah. 



(19) 



where each /C/ consists of products of such kinematic factors while the respective integrals 
can be written as 



nab 



= / dxi... 
Jo Jo 



1 Ar-3 
dXN-3 



n 

a=l 



N-3 2a'fc(,+3( 



a 



b=a 



2a' k2+ak3+b+nab 



(20) 



j=a 



with the indices b > a — 1,2, . . . , N—3, and the integers Ua, riab taking values 0, ± 1 or ±2. 
By convention, the sum in the exponent is zero for b — a. The integral involves N{N — 3)/2 
different Laurent polynomials in Xa- Their integer powers Ua, Uah control the physical poles 
of the amplitude, in N{N — 3)/2 invariant masses of dual resonance channels involving 



2, 3, ... , E{^) external particles {E denotes the integer part). 
For iV = 4, the integral (20) yields the Beta function 



m 
nil 



= f dXi xf (1 - Xi)^i^+'*" = — 

Jo 



•323+^1 — 1 , — S12— nil 
S23+ni 



S23 + ni - 1 



(21) 



For N — 5, one obtains the hypergeometric function 3F2 [17]: 



ni,n2 
nii,ni2,n22 



1 1 

J dxi J dx2 Xi 



3+S23+ni „-2+si5+n2 



X (1 - (1 - X2)'*^+"'' (1 - XiX2)^^^+"^' 

r(s23 + ni - 2) r(si5 + n2 - 1) r(s34 + nil + 1) r(s45 + «22 + 1) 



r(s23 + S34 + ni + nil - 1) r(si5 + -545 + n2 + 1122) 

X 3F2 



(22) 



S23+"l-2 , Sis+n2-l , -S35-ni2 . 
S23+S34+ni+nil— 1 , S15+S45+n2+Tl22 ' 



Both integrals (21) and (22) boil down to hypergeometric functions of one variable, i.e. 
some pFq = 1 . However, this pattern does not persist beyond = 5, due to 

the form of the integrand (20) that does not fit into any hypergeometric function of one 
variable u. In general, one obtains multiple Gaussian hypergeometric series, more precisely 
certain generalized Kampe de Feriet functions [18]. For example, for N — 6, the integral 

111 

F \ "^'"^'"^^ ] = [dxJ dX2 I dXs a;r'+^^^+"^ S+a'(kMk,)Hn2 -2+.i6+n3 

[nii,ni2,n22,ni3,n23,n33 J III 



X (1 - xi)^3*+''" (1 - X2)^*^+"'' (1 - a;3)^^«+"33 (23) 

X (1 - XiX2)'''+"^' (1 - X2X3y'''+'''' (1 - XiX2X3y+'''' 

can be expressed in terms of the triple hypergeometric function F*^^^ [19]. 

A very important part of the discussion of scattering amplitudes is the examination of 
their low-energy behavior. To that end, the integrals (20) must be expanded in powers of 
a'. One can first expand the integrand and then integrate the series term after term. A 
typical, but by far not the most general, class of integrals that appear in this way are: 

c(«i, fn ^-^rv-] f • • /' n ^^^^^ ■ (24) 

V,=i J Jo Jo 1-tlx. 

i=l 

They integrate to multiple zeta values of length k [20] : 

fe ^ 00 k / k 

c(.i,. ..,..)= Yi n E n E^o ' (25) 

ni>...>nfc>0 j=l 3 ni,...,nk=l j=l \i=j / 

with si > 2 , S2,...,Sk > 1. Such integer series are completely sufficient for discussing 
the expansions of amplitudes involving four and five gluons however, as mentioned before. 



starting ai N — 6, more general classes of integrals appear. Their expansions involve multiple 
harmonic series and generahzed Euler-Zagier sums. We refer interested readers to Ref.[19] 
for a detailed account on the relation between multiple Gaussian hyper geometric functions 
and Eulcr-Zagier sums. Actually, the integer sums that appear in the context of multi-gluon 
string scattering play an important role in modern number theory [21]. 

The number of independent (with respect to the momentum conservation constraint) 



kinematic factors ICj and of the associated functions F entering into the A'"-gluon par- 

L"a6 J 

tial amplitude (19) grows with N. In our analysis, we encountered 77 functions ior N — 5 
and 1, 270 functions for N — 6, although these numbers may vary depending on the imple- 
mentation of momentum conservation constraints etc. In fact, many functions are related 
by means of polynomial relations of their integrands or by partial integration. The only 
systematic way of handling them for arbitrary N is to find a basis, consisting of an a priori 
unknown number i/^ of functions, and to express all other functions as linear combinations 
of the basis elements with the coefficients given by some rational (homogeneous) functions 
of the kinematic invariants (5). This program has been successfully implemented in [17, 22] 
for N = 5 and in [19] for N = 6 and will be continued in [23]. For a given N, an efficient 
way of generating systems of equations relating the integrals (20), that can be used to find a 
minimal set of independent functions, is based on world-sheet supersymmetry [19]. It works 
in the following way. In Eq.(15), the two vertices in the (— l)-ghost picture were inserted, 
for convenience, at zi and Z2- However, due to world-sheet supersymmetry, they could be 
inserted at any other two points, hence there are (^) ways of computing the same amplitude 
that should give the same answer for the coefficients of all (independent) kinematic factors )Ci 
in Eq.(19). By comparing these coefficients, one obtains many relations among the integrals 
(20). The corresponding set of equations is always under-determined and may be solved by 
expressing all functions (20) in terms of a i/^v-dimensional basis. Of course, the dimension 
of the space of functions grows with the number of gluons: — 1, — 2, vq — 6, ... 

Although only one partial amphtude ^4(1^1, ... , N^'^) has been discussed here explicitly, 
all other partial amplitudes A{a{l^'^), . . . , a{N'^''^)) can be obtained in exactly the same 
way. As we shall see in the following, a convenient choice of i^n basis functions is dictated 
by various physical properties of the amplitude (1). 
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III. FIVE GLUONS 



The purpose of this section is to summarize the results of five-gluon computations [17, 22] 
and to rewrite the five-gluon amplitude in the four-dimensional helicity basis. Recall that 
up to five gluons, the amplitudes are purely MHV. Here, five invariants are necessary to 
specify the kinematics. They can be chosen as = [i]^, i = 1, . . . , 5, i.e. as the cyclic orbit 
of [Ijj^ obtained by the action of Z5 subgroup of cyclic permutations, generated hy i ^ i + 1 
mod 5 [24]. 

The integrals over two vertex positions have the form (22), specified by five integers 
ni,n2, nii,ni2, n22- One finds [19, 22] that all integrals can be expressed in terms of just two 
functions: 



/i 



2,1 
0,0,0 



and /: 



2 



3,2 
0,-1,0 



(26) 



By means of simple algebraic operations and partial integrations it is easy to see that these 
functions transform in the following way under the Z5 generator i ^ i + 1 mod 5: 

fi ^ F Q = ^— [ S2S5 fi + (S2S3 - S3S4 - S1S5 + S4S5) /2 ] , (27) 

L ' ' J S1S3 

/2^/2. (28) 



In the notation of [22], the partial amplitude 

A{1^\2^\ 5^5) = T-Aym{1^\2^\ 5^'^)+K3-Af4{1^\2^\ 3^\4^\ 5^^) (29) 
where 

T{si) = S2S5 fi + (S2S3 + S4S5) /2 and K-i{si) = /s . (30) 

In Eq.{29), Aym is the tree-level Yang- Mills amplitude while Ap4 is generated by the Xp4 
interaction term (12) discussed in the Introduction. There are two Feynman diagrams, 
shown in Figure 1, that combine to Ap4: the diagram with Ip4 four-gluon vertex including 
one off-shell gluon decaying into two external gluons via the standard three-gluon Yang-Mills 
interaction, and the diagram with Xp4 five-gluon vertex. The function in Eq.(29) can 
then be interpreted as a string "formfactor" of a'^Tr interactions, playing role similar to 
the Yang-Mills formfactor T. 

The amplitudes Aym and Api, as well as the functions T{si) and K^lsi) are invariant 
under cyclic permutations, therefore the amplitude (29) is cyclic invariant. Furthermore, 
it has the correct factorization properties into four-gluon amplitudes, in the limit of zero 
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FIG. 1: Feynman diagrams contributing to the Ap4 part of the five-gluon ampUtude involve a 
single four- or five-gluon vertex due to the 2p4 effective interaction, represented here by the blob. 



(soft) momentum of one gluon and in the limit of two parallel momenta [22]. The low-energy 
behavior of the amplitude is determined, up to the order 0{a'^), by the following expansions: 

Jl — 7- — \ — 

S2S5 6 \S2 S5 

+ C(3) -S1 + S3 + S4H \ \ \ +... , (31) 

V ^2 -55 S5 S2 J 

/2 = y - C(3) (si + S2 + S3 + S4 + S5) + . . . . (32) 

In order to rewrite the amplitude (29) in the MHV form, we evaluate it for the spe- 
cific configuration of the polarization vectors, choosing a gauge with the most convenient 
"reference momenta" [4, 5]. We choose the reference momenta for ^~(1,2) and ki for 
r(3,4,5): 

C(z) = _iMM£lforz = l,2 and tAj) = ^^^1^ for j = 3, 4, 5. (33) 

Indeed, with such a choice, the only non- vanishing scalar products of the polarization vectors 
are 

r(^)-^^<^)-i|{i| - ™.r(4).-MM. (34) 

In this gauge, the respective kinematic factors /C/, see Eq.(19), contain only one fac- 
tor while the remaining three polarization vectors are contracted with the momenta. The 
computation consists of manipulations with spinor products, involving a repeated use of 
the momentum conservation law and of Schouten identity [4, 5]. A very useful check is 
provided by the cancellation of unphysical poles [Si]"^ and (Ij introduced by the choice 
(33) of the reference momenta. After factorizing out the Yang-Mills MHV amplitude ^ym-> 
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c.f. Eq.(2), the remaining spinor products can be either expressed in terms of kinematic 
invariants Si or they form the products [4, 5] 

a'^T[i,j,l,m] = a''(ii)[jZ](Zm)M=Q;'2trQ(l-75)|(fi^,Jfi^^^ (35) 

= ^ [ % Sim - Sil Sjm + Sim Sji - 4i e{i,j, l,m)] . (36) 

The new feature, as compared to four gluons, is the appearance of Levi-Civita pseudoscalars. 
They originate from the Api part of the amphtude only. For five gluons, the momentum 
conservation law allows expressing all such pseudoscalars in terms of one of them, that can 
be chosen to be e(l, 2, 3, 4). The final result is 

A(l-,2-,3+,4+,5+) = [ V^'\s,)-2z P^'\s,) 6(1,2,3,4) ] OJlg,, (37) 

where 

y('5)(s.) = /i + i (S2S3 + S4S5 - S1S2 - S3S4 - S1S5) /2 and P^^\si)=f2. (38) 

The above functions, as well as the pseudoscalar e(l,2,3,4), are invariant under cyclic 
permutations [25], thus the factor multiplying the Yang-Mills amplitude in Eq.(37) is cyclic 
invariant. 

The low-energy behavior of the amplitude (37) is determined by the expansions (31): 

V^'\s,) = i-^{W 

+ ^ ( {sls2} + {sisl} + ) + ..., (39) 

P^'\si) = ^-C(3) {s^} + ... , 

where the curly brackets enclosing kinematic invariants imply the summation over all distinct 
elements of the respective cyclic orbit [26]. 

The connection to the four-gluon amphtude (9) can be estabhshed by considering the soft 
limit, say one ki — > 0, see Appendix A. Then the pseudoscalar part of the factor disappears 
due to the momentum conservation while the function 

^ h=o r(i + S1 + S2) ^ ' 

reproduces the Veneziano formfactor in Eq.(8). 
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All other partial amplitudes 74((t(1""), (7(2~), (t(3+), (t(4'*"), (7(5+)), which according to 
Eq.(l) are necessary for constructing the full MHV amplitude ^'^'''''(1", 2~, 3+, 4+, 5+) are 
obtained from A{1~, 2~, 3"*", 4+, 5"*") by simply applying the coset permutations a to the right 
hand side of Eq.(37). 

IV. SIX GLUONS 

The step from five to six gluons is highly non-trivial. Even in QCD, the original cal- 
culation [27] used some of the most advanced tools available at that time, like extended 
supersymmetry [14], a special choice of the color factor basis etc. In addition to algebraic 
complications due to large numbers of Wick contractions and of the associated kinematic 
factors, there is a new physics element appearing at the six-gluon level: the scattering ampli- 
tudes allow also some non-MHV helicity configurations. Furthermore, each kinematic factor 
brings an integral over three vertex positions. The new challenge is to find relations between 
more than one thousand of such integrals and to express them in a suitable basis. Before 
discussing this problem, we review the six-particle kinematics (see also Appendix A), which 
also exhibits some new features as compared to the five-particle case. 

A. Six-Particle Kinematics 

In five and more dimensions, the number of independent kinematic invariants in a six- 
particle scattering process can be counted by using the momentum conservation law. There 
are nine invariants that can be grouped into two irreducible representations of the Zg cyclic 
group generated by i — > i -|- 1 mod 6: 

= H„ i = 1, . . . , 6, tj = lil^, j = 1, 2, 3, (41) 

i.e. the Zg orbits of |1]^ and |1]2 [24]. In four dimensions, however, these variables are 
subject to a fifth-order polynomial constraint [28] that reduces the number of independent 
invariants from nine to eight. This is due to the trivial fact that in four dimensions, at 
most four momentum vectors can be linearly independent. Then the columns of the five by 
five Gram matrix s built of the elements Sjj, 1 < i,i < 5, cannot be linearly independent, 
therefore det s = 0. 
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The vanishing of the Gram determinant is closely related to the following identity involv- 
ing the metric tensor g^^^, and the Levi-Civita tensor eapjs' 

One can eliminate one four- momentum, say /ce, by using momentum conservation, and define 
the following pseudoscalars: 

ei = e(2,3,4,5) 62 = e(l, 3, 4, 5) £3 = e(l, 2, 4, 5) £4 = e(l, 2, 3, 5) £5 = e(l, 2, 3, 4) , 

(43) 

and further define the five-component vector e = (ei, —62, 63, —64, 65). Then the identity (42) 
implies 

^ -Sl e2 - (Si + S2- ti) 63 - (S2 + S^-ti- ^2) 64 - (S5 + Sq - ^2) 65 ^ 
Sl ei S2 63 + (S2 + S3 - ^2) 64 + (S3 + S6 - ^2 - ^3) ^5 

v = e-s^ -[si S2 - ^i) ei - ^2 ^2 - S3 £4 - (53 S4 - ^3) £5 = , 

(S2 + S5 - ti - t2) ei (S2 S3 - €2 + S3 €3 S4 €5 
\ -(S5 + S6 - ^2) ei - (S3 + S6 - ^2 - ^3) £2 - (S3 + S4 - ^3) €3 - S4 £4 y 

(44) 

where we introduced the vector v = {v^, v^, v^, v^, v^) with all components vanishing due 
to e • s = 0. Thus the vanishing of the Gram determinant, det s = 0, ensures self-consistency 
of the above identity. Although Eq.(44) will be important for understanding how six-gluon 
amplitudes transform under cyclic permutations, it is convenient to keep as many scalars and 
pseudoscalars as allowed by momentum conservation, without using the Gram determinant 
constraint or Eq.(44) explicitly to eliminate the redundant invariants. 



B. Integrals and Their Six-Element Basis 



The integrals (20) over three vertex positions have the form (23), with nine integers 
ni, ^2, 713, nil, ni2, ^22, 7^13, 7^23, ?T'33- Now six functions are necessary to form the integral 
basis. A convenient basis to start with is: 



Fi= F 
F2= F 



3,2,1 
0,0,0,0,0,0 



4,3,1 
0,-1,0,0,0,0 



F^= F 
F^= F 



4,3,2 
0,0,0,-1,0,0 



4,4,2 
0,-1,0,0,-1,0 



F,= F 
Fe= F 



4,3,2 
0,-1,0,-1,0,0 



4,3,2 
0,0,0,-1,-1,0 



(45) 
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In choosing the above functions, we were guided by their low-energy power expansions in a', 
by their soft hmits and by their transformation properties under cychc permutations. The 
a'-expansion of Fi, derived in Appendix B, is: 

_l_ ^^2^ / ^4 + ~ ^1 _|_ -53 + ^4 — ^3 _|_ -§4 + 34-2 ^ + 55^6 _|_ ^2S3_+^ j _^ 



\ S2 Sq S2SQ ^2^2 "56^2 / 

while the expansion of F2, see also Appendix B, starts with a single pole: 

f,.!^i_C(3)fi±ii±ii±^ + ... . (47) 

D Sq Sq 

These functions are related by the soft limit /cq — [29] to the five-gluon functions /i and 
/2of Eq.(26): 

5256^2-^1 > S2S5/1 , SeF2 > f2- (48) 

fc6=0 K6=0 

In fact, the expansions (46) and (47) are very similar to (31) and (32), respectively The 
remaining four functions have no poles. In particular, the low-energy expansion of F^ begins 
with the constant C(3): 

^3 = C(3) - ^C(4) {si + 4S2 + 3S3 + 2S4 + 3S5 + 4S6 + ti + 4t2 + t3) + • • • • (49) 

This function is not cyclic invariant; three additional functions, F4, F5 and Fg, are necessary 
in order to form a closed representation of Zg. Under the generator i ^ i + 1 mod 6, they 
transform in the following way: 

Fa^-Fg + Fs , Fe^Fs , F5^F4 , F4^F6. (50) 

The low-energy expansions of the functions F4, F5 and Fg also begin with C(3): 

F4 = 2C(3) - ^C(4) (7si + 5s2 + 5s3 + 7s4 + Ssg + 5sq + 5ti + 2t2 + Bt^) + ... , 

F5 = 2C(3) - ^C(4) (5si + 5s2 + 7s3 + 5s4 + Ssg + Tsg + 2ti + 5^2 + 5^) + . . . , (51) 

Fq = 2C(3) - ^C(4) (5si + 7s2 + 5s3 + 5s4 + Tss + 5sq + 5ti + 5^2 + 2^3) + . . . . 
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For completeness, we list here also the cychc transformations of Fi and F2. 

-S3S6 (S4 + S5- ti) (F4 - Fe) - Ssts {S2 + S5-ti- ^2) (-^3 - -^5) 

-56^2 {Si + S4-ti- ts) F3 - (siSe - S5S6 + 55^3 - t2t3) X (52) 

X [(S4 + S5- ti){F^ - Fe) + (si - S3 + S5 - ti){Fs - F5) + (s2 + S4 - - ^1)^3] , 
S6 F2 seF2 + S6 (F4 - Fe) + (s2 - h) (F3 + F4 - F5) 

+ (55 - t2) (F3 + F4 - F5 - Fe) . (53) 

Although all six-gluon integrals can be expressed in terms of the basis (45), we will see that 
the actual amplitude involves certain combinations that assume a simpler form when written 
in the original notation of Eq.(23). In Appendix C, we will express the relevant integrals in 
terms of F^, A; = 1, . . . , 6. 



C. MHV Amplitude 

The results of [19] allow expressing the full six-gluon string amplitude in terms of six gen- 
eralized hypergeometric functions (45), with each function multiplying a long combination 
of kinematic factors involving all possible contractions among the momentum and polariza- 
tion vectors. Unlike in the fivc-gluon case, the amplitude cannot be simply separated into 
parts associated to some functions like T and K3 that can be attributed to distinct effective 
interactions, see Eq.(29). Experience with QCD suggests that the complications are related 

to the existence of the non-MHV part, with ( 1" ++) hehcity configurations. It is 

reasonable, however, to expect that the ( h + + +) MHV amphtude can be simplified. 

To that end, we substitute to the general expression [19] the following polarization vectors: 

^-(i) ^ _(MM^ for i = 1, 2 and CU) = ^ for j = 3, 4, 5, 6. (54) 

The kinematic terms surviving in such a configuration contain only one ^i^j factor, and the 
remaining four polarization vectors are contracted with the momenta. A generic term has 
the form 

(12)(lj)(lA:)(lm)(ln) L7-6][fc6][a6][/jm][7n] ^^^^^^ „. ^ ^ ^ „, o 4 , . ^ „ 904, 

(13)(14)(15)(16) M26^ '^^'^ a^/3^7 = 3,4,5, j,k,m,n ^ 2,3,4,5 

(55) 
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times a linear combination of the six basis functions, with the coefficients being rational 
functions of scalar invariants. There are more than one thousand of such terms, so it is 
quite a tedious task to simplify the answer. The final result can be written as 

k=5 



^(l-,2-,3+,4+,5+,6+ 
with the functions: 



k=l 



(6) 

YM ' 



(56) 



p(6) 


= Si 


F 






+ 


p(6) 


= S2 


F 


p(6) 


= S3 


F 


p(6) 


= S4 


F 


p{6) 


= S5 


F 



4,3,2 
0,-1,0,-1,-1,0 



+ (S2 + S5-ti- t2) F 



4,3,2 
0,-1,0,-1,0,0 



3,3,2 
0,-1,0,0,-1,0 

4,3,2 

-1,0,0,0,-1,0 

4,3,2 
0,0,-1,-1,0,0 

4,3,2 
0,-1,0,0,0,-1 



4,4,3 
0,0,0,-1,-1,0 



+ (S3 + S6 - ^2 - ^3) F 
+ (Si + S4 - ti - ^3) F 
+ {S2 + S3 - S4 - ^2) F 

+ (S3 + S4 - S5 - ^3) F 



4,4,2 
0,-1,0,0,-1,0 

4,4,3 
0,0,0,-1,-1,0 

4,3,2 
0,0,0,-1,0,0 

4,3,2 
0,-1,0,-1,0,0 



(57) 



^ s^s^t2 F 



3,2,2 
0,0,0,0,0,-1 



+ 2 (*2S3 - S3S4 + S3S6 + 8^2 - S2t3 - t2h) 



(6) 



+ ^ (-S2S3 + S1S4 - S4S5 - S3S6 + Ssti - S4t2 + S2t3 + Ssts " ^1^3 + ^2^3) -^i^^ (58) 



+ 7: (S2S3 - S1S4 + S2S5 + S3S6 + S5S6 - ssti - seh - S2t3 - s^ta + tit2 + ^1^3 - ^2^3) -R 



(6) 



+ o ("■52S3 + S1S4 - S2S5 - S1S6 + S3ti + Soti + Sit2 + S2t3 " ^1^2 - ^1^3) 



(6) 



+ 2 {~^l^2 + S2S3 + S2S5 — Ss^i — Sit2 + ^1^2) P^^ — S5S3 P^^ + S5(S3 — ^2) -^3°'' ■ 

The result can be expressed in the basis of functions introduced in Section IV B by using 
the formulas written in Appendix C. 

Although the six-gluon V and P functions appear complicated, they have very simple 
transformation properties under cyclic permutations. After expressing them in terms of the 
basis functions F^, /c = 1, . . . ,6, see Appendix C, and using the transformation properties 
(50, 52, 53), it is easy to see that V'^^\si,ti) is cyclic invariant. Furthermore, the functions 
P^^^(si, tj) transform among themselves in such a way that the imaginary part of the Yang- 
Mills formfactor in Eq.(56) is also invariant. This can be seen in the following way. Let us 



,(6) 



j(6) 
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put the five functions into the vector P^^) ^ (p^(6)^ _pi6)^ ^ _pj6)^ p{6)y ^j^^^ 
the action of i — > i + 1 mod 6 on e and P^^^ can be written as 



e — ^ e M , P(^) — ^ P(^) (M*)-^ + A P3 , 



(59) 



with the unimodular matrix M and the vector A: 

/ 1 1 1 1 l\ 
-10 
0-10 
0-100 
\ -1 0/ 

However, with Eq.(44), we find: 



A* = 



f-Sl - S2 + S6 + 
-Si + S2 — Sq + 
Si + S4-ti- ts 
S3 - S4- S5+ ti 
\-S3 - S4 + S^+tsJ 



(60) 



€ A' = v' + + v^ + v'' = , 



(61) 



thus P^^^ e* P(^^ e*. As a result, similarly to the case of four and five gluons, the full 
string formfactor of the MHV six-gluon amplitude (56) is cyclic invariant. It also has the 
correct soft limits [29] when any momentum goes to zero: 



ki=0 



(62) 



J2i-^y^' Pi^'\s.U)-—^P^'\s,), P['\s,,U) P^'\s,) for l<5. 

'—^ fe6=0 fei=0 



l=\ 



Furthermore, it has the right collinear limits [29], when the momenta of adjoining gluons, 
ki and /cj+i, with i + 1 mod 6, become parallel: 



J2^kPlf\s,,U) — e(l,2,3,4)p(5)(s, 
fc=i 



(63) 



The low-energy behavior of the amplitude (56) is determined, up to the order 0{a ), by 
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the following expansions: 

y {Si, U) ^ l-'^i {S,S2} - {S,S,} + {hh} ) + ^ {{S,sl} + {SIS2} - {sis,} 

+{siS2^i} - {siS4^i} - {s2S5ti} " 3{siS4i2} + {si^i^s} + {titj} + 0?t2} + 3 ti^a^s) , 

2 

Pi'^\si, ti) ^ ^ + C(3) (Si + 2S2 -S3-S4 + 2S5 + S6- 3h - 3t2 - ts) , 
2 

{si, ti) ^ ^ + C(3) (2s2 + 2s3 - S4 - S5 + S6 - ^1 - 3^2 - 2^3) , (64) 

2 

Pt\si, ^0 ~ Y + C(3) (2S3 + S4 - S5 - S6 - - ^2 - 2^3) , 

^i^H^i, ^i) ~ Y + C(3) (-si + Ss + S4-Se-ti-t2- ts) , 
2 

Pi''\si, ii) ~ Y + C(3) {-si - S2 + S3 + 2s4 -ti-t2- 2ts) . 

All other partial amplitudes A{a{l'~), cr(2~), (7(3+), (j(4+), (7(5''"), (7(6+)) can be obtained 
from 74(1", 2~, 3+, 4+, 5+, 6+) by simply applying the S^j^L^ coset permutations a to the 
right hand side of Eq.(56). 



D. Reconstructing the Amplitude from First Principles 

After checking that the amplitude (56) satisfies all self-consistency conditions following 
from cyclic symmetry and soft/coUinear limits, we would like to proceed in reverse, in order 
to understand to what extent the form of the string factor is determined by these conditions. 
To that end, we make the following ansatz for the function y(^)(sj,ii): 

6 

y(^) (Si, ti) = S2 Se t2 Fi + (S6 F2) + Yl ' (^5) 

l<i<j<9 1=3 l<i<j<k<9 

in the basis of six functions (45), with 45 and 660 real constant coefficients Xij and Xlj/,, 
respectively. We will try to fix these constants by imposing the above self-consistency 
conditions. For convenience, we use the notation 57 = ^1, sg = t2, Sg = t^. The leading 
term in the a' expansion of the ansatz (65) which, according to Eq.(46), is equal to 1, is 
dictated by the zero slope (Yang-Mills) hmit, V'^^^Si^ti) ^ 1. The next-to-leading order 
C(q;'^), with the common C(2) = ^ factor, is governed by the constants Xij. They are 
completely determined by the cyclic symmetry and soft limits. Finally, after imposing the 
right coUinear limits, all remaining constants j^. can be expressed in terms of one of them. 
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Aygg. In this way, the ansatz (65) becomes 

V^''\si,U) = V^'\s„U) + A^89 ( ht2ts - { sis^t2 } ) (F4 + Fs + Fe) , (66) 

thus the real part of the string factor is completely determined by the self-consistency con- 
ditions, up to one constant. The respective term is cyclic invariant and vanishes in both soft 
and coUinear limits. 

In order to examine the imaginary part of the string factor, we assume that it has the 
form emPm {si,tj), with the following ansatz 

m=l 

6 9 

Pi^) {S„ U) ^SeF^ + J^FlY. ^''rnk Sk (67) 

1=3 k=l 

for the five functions Pm\si, ti). Here again, we try to fix 180 real constants by demand- 
ing that the sum ^ Pm be cyclic invariant and that it has the correct soft/coUinear 

m=l 

limits. The latter requirement fixes 156 of 180 constants. By further imposing cyclic invari- 
ance, one ends up with only four arbitrary constants, //|g, /if^, /if^ and //fg. Finally, after 

using the relations between pseudoscalar invariants, written as f = in Eq.(44), one finds 

5 5 

m=l m=l 

- (3 4^ + 1) [{se - t2) €5 + (S2 - ^2) £4 + S2 63] {Fe + F5 - F4 - 2 F3) , 

thus also the imaginary part of the string factor is determined up to one constant. 

There is one more constraint available. In the Abelian case, the leading term in the a' 
expansion must vanish because it is entirely due to Yang-Mills gluon self-interactions. The 
next-to-leading term, which is associated to the Xp4 interaction, must also vanish: while 
four Abelian gauge bosons interact via the corresponding Born-Infeld term, they cannot 
spread via Yang-Mills interactions, like in the left diagram on Figure 1. By requiring that 
the a'-expansion of the Abelian amplitude starts at order higher than 0{a''^), we obtain 
A^gg = and //fg = — 1/3 [23]. Actually, with this constraint, the Abehan amplitude starts 
at order 0{a"^) [with the common factor of C(4)]- We conclude that the six-gluon MHV 
amplitude can be uniquely determined from first principles. It is worth mentioning that 
in the simpler five-gluon case, the cyclic symmetry and soft limit are completely sufficient 
to determine the amplitude, similarly to the constants Ajj that govern the next-to-leading 
0{a''^) contribution to the six-gluon amphtude. The reason, to be elaborated in the next 
section, is that all 0{a''^) terms originate from 2^4 interactions. 
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V. N GLUONS 



It is clear from the discussion oi N — 5, and especially oi N — 6, that the computational 
complexity increases steeply with N. The integrals (20) become more comphcated and the 
number of independent functions grows. The functions emerging in the step from A^— 1 
to N have low-energy expansions starting at 0{a'^~^) [with a common factor of ('(A'^— 3)]. 
On the other hand, the simple, factorized forms of Eqs.(9), (37) and (56) strongly suggest 
that MHV configurations enjoy a special status. The fact that = 5 as well as N — 6 
MHV amplitudes can be reconstructed from first principles, by using very simple physical 
constraints, is very encouraging because it opens way to an iterative procedure suitable for 
larger numbers of gluons. It also indicates the existence of some recursion relations similar 
to those in QCD [7]. A recursive construction of the amplitudes requires however a better 
understanding of the space of generalized hypergeometric integrals (20). This ingredient 
will have to wait until completion of Ref. [23]. Nevertheless, already at this point, we can 
determine the leading 0{a''^) string corrections. 

An A'"-gluon scattering process can be parameterized in terms of A'^(A'" — 3)/2 kinematic 
invariants which can be chosen as the cyclic orbits of k — 1, . . . ,-E(y — 1), where 

E denotes the integer part. Recall that the cyclic Zjv group is generated by the shift of 
indices labeling gluons from i ^ i + 1 mod A^. Note that for A^ odd, the last orbit contains 
A^ elements, while for A^ even their number is reduced by the momentum conservation to 
A^/2. As in the case oi N — 6, we can ignore the four-dimensional Gram determinant 
constraints [28] that reduce the number of independent invariants to 3A'" — 10. We also keep 
the pseudoscalars e(A;, Z,m, n), with k<l<m<n<N, which are independent as far as 
the momentum conservation is concerned but are related by equations similar to (44). 

Let us first collect all known 0{a''^) terms and rewrite the leading and next-to- leading 
terms in the low-energy expansions of A'" = 4, 5, 6 MHV amplitudes as 

' fmi.^ + 0{a'^) , (68) 



12 ^ 



/l(l-,2-,3+,4+,---,Ar+) = 
where Q^^^ are the following Lorentz- invariant, homogenous of degree four, functions of the 
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FIG. 2: Feynman diagrams contributing the leading 0{a''^) string corrections to A^-gluon YM 
amplitudes involve one Ipi effective interaction vertex, while the remaining vertices are due to the 
tree- level YM interactions. 

momenta: 

Q(^) = SiS2 + S2S3 + S3S4 + S4S5 + S5Si+4ze(l,2,3,4) (69) 
Q('5) = S1S2 + S2S3 + S3S4 + S4S5 + S5S6 + SqSi + tit2 + t2h + ^3^1 - SlSi - S2S5 - S3S6 

+ 4t [e(l, 2, 3, 4) + e(l, 2, 3, 5) + e(l, 2, 4, 5) + e(l, 3, 4, 5) + e(2, 3, 4, 5)] 

At this order, the corrections are generated by gluonic tree diagrams involving only one Xp4 
vertex, c.f. Eq.(12), and a number of standard Yang- Mills interactions. A typical diagram 
contributing to A^-gluon scattering is shown in Figure 2. In order to determine Q'-^-* for 
arbitrary one can either calculate the sum of such Feynman diagrams or apply iteration 
utilizing the soft limit and Z^r symmetry. Although it is possible to formulate the latter as a 
formal recursion relation, we prefer to apply the iterative procedure explicitly, step by step, 
starting from Q^'^\ Here again, it will be very convenient to use the notation introduced 
in Section III: an expression enclosed inside curly brackets {■ ■ ■ } denotes the sum over all 
distinct elements of its cyclic permutation group orbit. Thus, for example, 

Q^'^^ = {S1S2} + {tit2} - {S1S4} + J2 <k,l,m,n). (70) 

k<l<m<n<6 

For N = 7, the 14 invariants are = [ijj^, i = 1, . . . ,7 and tj = [jjg, j = 1, . . . , 7. The 
real part of Q*-^-* must be a cyclic invariant, quadratic form in Sj, tj. There are 15 quadratic 



23 



cyclic invariants, but there exists only one linear combination, 

ReQ(7) = {siS2} + {tit2} - {sih} (71) 
that gives ReQ^^^ in the soft limit 

fcr ^ : sq^O, sr ^ 0, ^ h, h ^ s^, U ^ sq, tr ^ si. (72) 

In order to determine the imaginary part, we first note that although there are 15 linearly 

independent pscudoscalars, there are only 3 Zy-invariant combinations. We list them below, 
together with their /cy — > limits: 

e(l, 2, 3, 4) + e(l, 2, 3, 6) + e(l, 2, 5, 6) + e(l, 4, 5, 6) + e(3, 4, 5, 6) ^ e(2, 3, 4, 5), 

e(l, 2, 3, 5) + e(l, 2, 4, 5) + e(l, 3, 5, 6) + e(2, 3, 4, 5) + e(2, 4, 5, 6) e(l, 3, 4, 5) + e(2, 3, 4, 5), 

e(l, 2, 4, 5) + e(l, 3, 4, 5) + e(l, 3, 4, 6) + e(2, 3, 4, 6) + e(2, 3, 5, 6) e(l, 2, 3, 5) + e(l, 2, 4, 5). 

Thus the imaginary part of Q^^-* is also uniquely determined to be: 

ImQ(^)=4 J2 (^{k,l,m,n). (73) 

k<l<m<n<7 

The above iteration can be continued to a larger number of gluons, with the unique answer: 
Q^^^ = 1: {[11J21J- fl {[llj2l,_J + CW+4i 6(fc,Z,m,n), 

fe=l fe=3 k<l<m<n<N 

(74) 

-{Ilk-2lf + lk_2} A^>4,even 



-{[ll-=5[^]^} TV > 5, Odd. 



It is very interesting that Eq.(74) bears a striking resemblance to the one-loop all positive 
hehcity amplitudes of QCD [30] . The resemblance originates at the five-gluon level, c.f. tt^ 
terms in our Eq.(39) vis-a-vis Eq.(5) of Ref.[30], and then propagates to N gluons because 
in both cases the multi-gluon results are uniquely determined by the permutation symmetry 
and soft limits. 

Work on the recursive construction of MHV amplitudes to all orders in a' is in progress 
[23]. 



24 



VI. CONCLUSIONS 



The main result of this paper, in addition to the specific formulas for scattering ampli- 
tudes, is the demonstration of a striking simplicity hidden in multi-gluon scattering, per- 
sisting at the full-fledged string level. The maximally helicity violating configurations retain 
their special status even after the tree diagrams are replaced by a disk world-sheet. The 
string effects are succinctly summarized in a number of kinematic functions, extending the 
well-known result for four gluons to an arbitrary number of gluons. We argued that the soft 
and coUinear factorization properties, combined with the Abelian limit, are completely suf- 
ficient to determine all A'"-gluon MHV amplitudes, however a completely recursive construc- 
tion requires a better understanding of the boundary integrals determining the kinematic 
functions. Work in this direction is in progress [23] . 

In superstring theory, it is often possible to describe a single physical process in several 
ways, by using various dualities. We believe that the simplicity of MHV amplitudes reflects 
the existence of an underlying integrable structure not only in QCD, but also in type I string 
theory. Hence it would be very interesting to understand if there is any room in the twistor 
formulation of string theory [8] that would allow accommodating open string corrections to 
YM scattering amplitudes. 

Another duality relevant to the present work is the type I-heterotic duality [31, 32]. Here, 
the disk-level interactions of 2n gauge field strength tensors in type I theory are dual to the 
heterotic {n — l)-loop interactions. In this context, it would be interesting to investigate 
a possible relation of our results to Ref.[30] and to the recent computations of all one- loop 
MHV amphtudes in QCD [33] . Our results should also help in explaining why the heterotic 
six-gluon amplitudes are not compatible at the two-loop level with a semi-classical type I 
description in terms of a naive extension of Born-Infeld electrodynamics [34, 35]. With 
some more work, the six-gluon amplitudes presented here could be used to extract the 
complete type I effective action and could shed more light on the long-standing problem 
how to construct a non-Abelian generalization of the Born-Infeld Lagrangian. 

When looking at the simple formulas describing multi-gluon superstring scattering, in- 
corporating infinitely many interactions among infinite number of particles, one cannot stop 
wondering if the effective field theory is really the right framework for describing low-energy 
string physics. For instance, the five-gluon amplitude (37) was used in [22] to determine 
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the complete 0{a'^) string corrections. The effective action consists of hundreds of terms 
and does not give justice to Eq.(37). The advantage of using the effective field theoretical 
description is that, in principle, it allows going off-shell and studying the modifications of 
classical field equations. However, it is certainly not the most efficient way of recording 
the S-matrix. Historically, string theory grew out from S-matrix theory but its formalism 
has evolved more and more towards Lagrangian quantum field theory. We need a better 
formalism, somewhere halfway between S-matrix and Lagrangians. 
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APPENDIX A: KINEMATIC INVARIANTS, SOFT LIMITS AND 
COLLINEAR LIMITS 

1. Kinematic Invariants 

The tables below contain the scalar products Sij = 2a' kikj, with i and j labeling rows and 
columns, expressed in terms of the kinematic invariants of type (5), used in the paper to 
describe multi-gluon scattering processes, for N — A, 5 and 6 gluons. 

N = 4 
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4 




1 


Sl 


-Sl - S2 




1 


2 






-Sl - S2 


2 


3 






3 



iV = 5 





2 


3 


1 


5 




1 


Sl 


— Sl — S2 + S4 


S2 - S4 - S5 


S5 


1 


2 




S2 


-S2 - S3 + S5 


-Sl + S3 - S5 


2 


3 




S3 


Sl - S3 - S4 


3 


4 




S4 


4 



iV = 6 





2 


3 


1 


5 


6 




1 


Sl 


-Sl — S2 + ti 


S2 + S5 - ti - t2 


-S5 — S6 + ^2 


S6 


1 


2 




S2 


-S2 - S3 + t2 


S3 + S6 - t2 - ^3 


-Sl - S6 + t3 


2 


3 




S3 


-S3 - S4 + ^3 


Sl + S4 - ti - ^3 


3 


4 




S4 


-S4 - S5 + ti 


4 


5 




S5 


5 



2. Soft Limits 



The soft limit is defined as /cj — > for some i. For a cyclic invariant function of the momenta, 
it is sufficient to consider only one soft momentum, say kN — > 0. Then the scalar invariants 
describing iV-gluon kinematics have the following limits in terms of the invariants describing 
A'^ — 1 gluons, for N = 5 and 6: 



N = 5 
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S2 


S3 


S4 


S5 






k^^O 


Sl 


S2 


Sl 








AT = 6 






Sl 


S2 


S3 


S4 


S5 


S6 




t2 


^3 


ke^O 


Sl 


S2 


S3 


S4 








S4 


S5 


Sl 



As k^ — > 0, the five-gluon Levi-Civita pseudoscalar invariant e(l,2,3,4) = £ — > 0. For 
N — 6, the soft hmits of pseudoscalar invariants defined in Eq. (43) are written below: 
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N = 6 





ei 


£2 


£3 


£4 


£5 












e 



3. CoUinear Limits 

The coUinear limit is defined as two adjoining momenta ki and /cj+i, with i + 1 mod A^, 
becoming parallel. Due to cyclic symmetry, these can be chosen as ^at-i and kN, with /cat-i 
carrying the fraction x of the combined momentum k^^i + ^ /^iv-i- Formally, 

kN-i xkN-i , kN ^ (I- x)kN-i , 

where the momenta appearing in the limits describe the scattering of A^" — 1 gluons. For 
N — 6, the coUinear limits of scalar invariants, written in terms of the invariants describing 
N — 5 scattering, are: 



N = 6 





Sl 


S2 


S3 


S4 


S5 


S6 


tl 


t2 


^3 


/cs — > X k^, /ce — > (1 — x)k5 


Sl 


S2 


S3 


XS4 





(1 - X)S5 


S4 


S5 


xsi + (1 - x)ss 



The coUinear limits of pseudoscalar invariants are written below: 



N ^6 





ei 


£2 


£3 


£4 


£5 


k^ ^ X k^, /cg — > (1 — x)k5 


e 


—xe 


xe 


—xe 


xe 



APPENDIX B: a' EXPANSIONS OF TRIPLE HYPERGEOMETRIC 
FUNCTIONS 

Most of the a' expansions of triple hypergeometric functions (23) presented in Ref.[19] 
apply to non-singular functions without poles, like F3, see Eq.(49). In that case, the ex- 
pansions of the integrals (23) can be directly mapped to convergent Euler-Zagier sums. In 
this appendix we derive the expansions (46) and (47) for the singular functions Fi and F2, 
respectively. 
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1. Fit Triple Hypergeometric Function with a Triple Pole 



Let us divide the integral defining Fi into two parts, Ii and l2'- 



3,2,1 
0,0,0,0,0,0 




{l-xy {l-yY* {l-zy 
X (1 - xyf^-'^-'^ (1 - yzf'-'^-'^ (1 - xyzy'+'*-^'-^' ^ + , (Bl) 



with: 



7i = (^J^dx x''-^ (1 - xY'^ (^j^dyj^dz y^'-^ z''-^ (1 - y)'^ (1 - ^)'« (1 - yzf'-'^-''^ , 



V 

. r(s2) r(i+s3) 

r(l+S2+S3) 



■ r(ii6) r(t2) r(i+«4) r(i+s5) 

r(l+S4+t2) r(l+S5+S6) 



Ht2' J 



/•I /.I /•! 

72 = dx dy dz x''-^ y^'-^ z''-^ (1 - x)'' (1 - y)^* (1 - z)'^ (1 - yzY'-'^-'' 
Jo Jo Jo 

X [ (1 - xyf^-'^-'^ (1 - - 1 ] ■ (B2) 



The first integral /i involves the Beta function (21) and the hypergeometric function 3F2 
(22). In fact, the latter integral is /i of Eq.(26) with appropriate arguments and its a' 
expansion can be found in Eq.(31). On the other hand, expanding the Beta function is 

straightforward, so altogether we obtain: 

r(s2)r(i + s3) 1 



r(l + S2 + S3) S2 



-C(2) S3 + C(3) S3 (S2 + S3) + ... 



(B3) 



r(s6) r(t2) r(i + s4) r(i + s5) 

r(l + S4 + t2) r(l+S5 + S6) 



3^2 



•56,^2, •S4+S5-*! . ^ 
l+S5+Se,l+S4+t2 ' 



^6^2 ^ Oe ^2 



+C(3) [s, + s,-t, + '' is. + h) ^s, (sa + sa) 



Se ^2 
Hence, we obtain the following a'-expansion for the integral Ii. 



+ ... . 



1 



52^6^2 



C(2) 



Si ^ _S5_ _^ S3 



S2S6 52^2 S6t2 



+ C(3) 



S4 + S5 - ^1 

S2 



+ 



Sdt2 , -4 + =^556 , S2S3 + S3 



S2S6 



52^2 



+ 



^6^2 



+ ... 



(B4) 
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The second integral I2 has a single pole in sq originating from 2; — > in the integrand. Its 
expansion in a' amounts to expanding it in powers of sq: 

si' : -(.3 + .4 - h) [dxfdy [dz .--^ = {ss + s,-ts) 

Jo Jo Jo S6 

_ S3 + S4-h f^^f^y H'^ - xy) -x) + S2 Ux) + S4 ln(l -y)+t2 Uy)] 

s& Jo Jo xy 

, (^3 + ^4- ^3)^ f\ f\ Hl-xyf C(4) (^3 + ^4- ^3) , ^ ^ ^,^^,1 

H 7, / dx / dy = 4 {S2 + S3 + S4 + ^2) + ^3 , 

2 S6 Jo Jq xy 4 S6 

: {si + S4-h- t^) fdxldy f dz ^"(^ ~ ^^^) ^ j^s^-h- t^) C(4) . (B5) 

Jo Jo Jo xyz 

Here, we have applied the following basic Euler integrals of the type (24): 

/■\ /"^ Inxln(l-xy) ,,,, [\ f\ ln(l - x) ln(l - xy) 5,,,, 

k/d, ^ i C(4) . = -C(3) . 



Jo xy 2 Jq Jq xy 

Cdxfdyfdz Mi:i^ = -C(4). (B6) 
Jo Jo Jo xyz 

In this way, wc obtain 

h = + " ''^ C(3)-^ + " [4 {s2+ss+s,+t2)+ts]-{s,+s,-t,-ts) C(4)+. . . 

(B7) 

Finally, after putting together (B4) and (B7), we obtain Eq.(46). 

2. F2: Triple Hypergeometric Function with a Single Pole 

Here again, we divide the integral defining F2 into two parts, Ii and l2- 
"1 i>i pi 

dx dy dz x'^ y*^ ^"^"^ (1 - xY' (1 - y)'' (1 - 



4,3,1 
0,-1,0,0,0,0 



'0 JO JO 

X (1 - xyY'-''-'*-' (1 - yzY'-''-'' (1 - xyzY'^"'~^'~*' = h + h (B8) 



with: 



h = (^j^dx z'^-' (1 - zY'^ (^jdxjdy x'^ (1 - xY^ (1 - yY^ (1 - xy)*3-^3-^*-^^ , 



-«5) rc^e) 

r(l+S5+S6) 



. r(l+sg) r(s6) _ r(l+32) r(l+s3) r(l+a4) r(l+t2) „ [l+sj, l+t2, l+S3+S4-*3 .il 

r(2+S2+S3) r(2+S4+t2) ^ [ 2+S2+S3, 2+S4+t2 ' J 



j-l j-1 j-l 

h ^ dx dy dz x'^ y*^ z'^'^ (1 - x)'^ (1 - y)'* (1 - zY^ (1 - a;?/)*^-'^-'*-^ 



10 Jo Jo 



X [ (1 - yz)*i-^^-^5 (1 - xyz)"i+"*-*i-*3 - 1 ] ■ (B9) 
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The first integral h involves the Beta function (21) and the hypergeometric function 3F2 
(22). In fact, the latter integral is /2 of Eq.(26) with appropriate arguments and its a' 
expansion can be found in Eq.(31). On the other hand, expanding the Beta function is 

straightforward, so altogether we obtain: 

^rnl^^l = - - C(2) .5 + C(3) .5 (.5 + ^6) + ... , (BIO) 

r(i + S2) r(i + S3) r(i + ^4) r(i + 12 



r(2 + S2 + S3) r(2 + S4 + ^2) 



■3-^2 



1+.S2, l+*2, l+S3+S4-t3 . 
2+S2+S3, 2+S4+t2 



C(2) - C(3) {s2 + S3 + S4 + t2 + h) + ... 



Hence, we obtain the following a'-expansion for the integral Ii. 

C(2) _ ^^3^ £2+^3+^4 + ^2 + ^3 ^ _ ^3^^^ 

The integrand of the second integral I2 remains finite for z — > 0. Up to the first leading 
order, it involves the following finite sub-integrals: 



/ Jo ''Jo ~ z{\-xy) 4 ' Jq Jq ''Jo " {1 - xy) z 4 

With this information we obtain 



{S4 + S5 - h) C(4) -^{si + s^-h- h) C(4) + . . . . (B12) 

After putting together (Bll) and (B12) we obtain Eq.(47). 

Finally, the four functions F3, F.^, F5 and do not contain any poles in the kinematic 
invariants (5). Hence their a' expansions can be obtained by the methods described in [19], 
i.e. by evaluating the relevant Euler-Zagier sums. 

APPENDIX C: FUNCTIONS V^^^ AND pj;^^ EXPRESSED IN THE BASIS 
Fk, fe = l,...,6 

The functions P^^^ and V^^'' governing the six-gluon MHV amplitude (56) are expressed 
in Eqs.(57) and (58) in terms of certain generalized hypergeometric integrals, in the notation 
of Eq.(23). In Section IV B, we introduced a basis of six functions, see Eq.(45), which is very 
convenient for studying cyclic properties and low-energy limits. The integrals appearing in 
Eqs.(57) and (58) can be expressed in this basis by using the relations obtained in Ref.[19] 
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as a combined result of partial integrations, use of world-sheet supersymmetry etc. The 
functions that enter Eq.(57) are: 



4,4,3 
0,0,0,-1,-1,0 

4,3,2 
0,-1,0,-1,-1,0 



S2 F 



S3 F 



S4 F 



S5 F 



3,3,2 
0,-1,0,0,-1,0 



4,3,2 
-1,0,0,0,-1,0 



4,3,2 
0,0,-1,-1,0,0 



4,3,2 
0,-1,0,0,0,-1 



- —-^3 + -^6 , 

= S6 (F2 + F4) + (.52 + S5 - tl - t2) (F3 + F4 - F5) 

+ (si - S5 - S6 + h) Fq , 

= S6 (F2 - Fs) + (s2 -S3 + S5-ti+ ta) (F3 + F4 - Fg) 

- {Sl -S3 + S5- ti) F5 + (Si + S2- ti) Fq , 

= S6 {F2 - F3) + {Si + S2- h) F3 - S3 {F3 - Fe) 

- (si -S3 + S5- h) F5 + (S4 + S5- h) (F3 + F4- Fe) , 
= se F2 + (s4 - S5- se + t2) F3 + (s4 + S5 - h) F4 

+ {si -S3 + S5- h) {F3 - F5) , 

= S6F2 + (Sl + S4-ti- t3) (F3 - F5) + (S4 + S5- h) F4 



The additional function that enters Eq. (58) is: 



S2S5 F 



3,2,2 
0,0,0,0,0,-1 



= S2S6 Fi - S2 (Sl - S5 - ^3) -^3 

+ (S4 + S5 - il) [ Se {F2 - F3) - {S3 -S5 + ti- t3) {F3 + F4) 
- {Sl - S3 + S5- ti) F5 + (Si + S3- S5- ^3) Fq ] . 
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